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STURMIAN MULTIPLE ZEROS FOR STOKES AND NAVIER-STOKES 
EQUATIONS IN M 3 VIA SOLENOIDAL HERMITE POLYNOMIALS 

V.A. GALAKTIONOV 

Abstract. The Cauchy problem for the 3D Stokes and Navier-Stokes equations, 
u t = -Vp + Au, divu = in M 3 x(-1,0], and 
,— i I u t + (u- V)u = -Vp + Au, divu = in M 3 x(-1,0], 

where u = [u, v, w] T is the vector field and p is the pressure, is considered. A smooth 
l/"~) ' bounded initial data Uo(a;), with divuo = 0, are prescribed at t = — 1. 

The problem of formation of multiple zeros at a given point (x,t) — (0, — ) of the 
components of u(x,t) is considered. In recent years, such a classic problem, which, for 
Ph ' the ID heat equation, was solved by Sturm in 1836, was under scrutiny for a num- 

ber of parabolic, hyperbolic, elliptic, and dispersion PDEs. As usual, such an analysis 
gives insight into a "microscopic blow-up scale properties" of the equations (j0.1[) under 
consideration. It is shown that formation of multiple zeros of solutions can follow "self- 
focusing" of nodal sets moving according to zero surfaces of the corresponding solenoidal 
Hermite polynomials as eigenfunctions of a rescaled adjoint Hermite operator. This is 
shown to be always the case for the first problem in (|0.1[) . but not always for the second 
one. Using such blow-up asymptotics allows to state a unique continuation theorem. 
lO \ A similar phenomenon is studied for the well-posed Stokes and Burnett equations with 

the bi-Laplacian — A 2 u instead of the standard viscosity operator Au in (|0.1|l . 

en 

1. Introduction: towards micro-scale structure of smooth solutions 

1.1. Stokes, Navier— Stokes, and Burnett equations. We consider the Cauchy prob- 
lem for the three-dimensional linear Stokes and Navier-Stokes equations (NSEs) 

(1.1) u t = -Vp + Au, divu = in R 3 x(-1,0], and 

(1.2) u t + (u- V)u= -Vp + Au, divu = in R 3 x(-1,0]. 

where u = [u,v,w] T is the vector field and p is the corresponding pressure. A smooth 
bounded initial data u (x) are prescribed at t — — 1, with divu = 0. 

The problem of formation of multiple zeros at the point (x,t) = (O.CD) of the vector 
field u(x, t) is considered. In recent years, such a classic problem, which, for the ID heat 
equation, was solved by C. Sturm (1836), was under scrutiny for a number of parabolic, 
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hyperbolic, elliptic, and dispersion PDEs; see short surveys for each type of linear and 
nonlinear PDEs in [17], and also [20], containing a most recent survey. As usual, such 
analysis gives insight into a "microscopic scale properties" of the equations ( 10. li under 
consideration. Indeed, such a microscopic blow-up approach was and is key in classic 
regularity problems including fundamental questions of the regularity of characteristic 
boundary points and singular points in potential and other related PDE theory; we refer 
to recent surveys in [T7J [19j [20] , where matching blow-up techniques were principally used. 

In the present paper, in a similar manner, we show that formation of multiple zeros 
of solutions can follow "self-focusing" of nodal sets moving according to zero surfaces of 
the corresponding solenoidal Hermite polynomials as eigenfunctions of a rescaled adjoint 
Hermite operator. It turns out that this is always the case for the Stokes problem. For 
the NSEs, the situation is shown to be more complicated, though similar zeros do exist. 
This allows us to state a non-standard unique continuation theorem for such problems. 

For more clear expressing our main "blow-up" techniques and their applicability in 
general PDE theory, we develop in Appendix A (the B one contains the corresponding 
Hermitian spectral analysis) at the paper end, as a natural extension, a similar multiple 
zeros analysis of the Cauchy problem for the fourth- order Stokes-like equations and well- 
posed Burnett equations: 

u t = -Vp-A 2 u, divu = in R Ar x(-l,0], u(x, -1) = u (x), 

u t + (u • V)u = -Vp - A 2 u, divu = in E Jv x(-l,0], u(x, -1) = u (x), 

where initial data Uo are sufficiently smooth and satisfy divuo = 0. Here, we have the bi- 
harmonic diffusion operators on the right-hand side of the u-equations. It turns out that 
our general scheme describing multiple zeros analysis can be applied, however, requiring 
another non-self-adjoint spectral theory for the corresponding rescaled operator, where 
generalized solenoidal Hermite polynomials naturally occur (Appendix B). 

1.2. Leray blow-up rescaled variables: why Hermite polynomials occur on 
micro-scales. For semilinear NSEs (II. 2p . we perform Leray's type \26y\ nonstationary 
blow-up scaling with the blow-up time T = 0: 

(1.4) u(x,t) = ^ t u(y,T), p(x,t) = j±-jP(y,T) } y = ^=, r = -ln(-f) -> +oo 

as t — y 0~. This yields the rescaled equations for u = (it 1 , u 2 , u 3 ) T and P, 

(1.5) u T = Au- \(y- V)u- |u- (u ■ V)u - VP, divu = in R 3 xR + . 

As a standard next step, we exclude the pressure from the equations (jl.5p . 



1.6) 



u r = H(u) = (B*- i/)u-P(u- V)u in 
where Pv = v - VA _1 (V • v) (||P|| = 1) 



P 3 



x 



In particular, Leray proposed not only to look at a self-similar blow-up as t — > T~ but also at a further 
similarity extension for t > T, i.e., in the complementary limit t —¥ T + , so that the blow-up factor \/T — t 
is replaced by \/t — T. See some historical and further comments on Leray's blow-up scenario of 1934 for 
the 3D NSEs can be found in [H § 2.2]. 
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is the Leray-Hopf projector of (L 2 (IR 3 )) 3 onto the subspace {w G (L 2 ) 3 : divw = 0} of 
solenoidal vector fieldsj. Another representation is Pv = (vi — R\a, v 2 — R20, V3 — R 3 a) T , 
where Rj are the Riesz transforms, with symbols Ifr, and a = R\Vi + R 2 v 2 + R3V3. We 

then first apply P to the original velocity equation in (II .2j) and next use the blow-up 
rescaling (11.41) . Using the fundamental solution of A in M 3 

(!-7) HV) = -ilJ[. 

the operator in (11.61) is written in the form of Leray's formulation [2H1 P- 32] 

H(u) = (B* - i l)u - (u • V)u + C 3 / r^Efs tr(Vu(z, r)f dz, 

(1.8) ' . R3 

where tr(Vu(z, r)) 2 = ^ (ij . } uj^ and C 3 = ^. 

It follows from (jl.8p that, in order to describe asymptotic behaviour of small solutions 
near multiple zeros, as a first step, one needs a spectral theory for the linearized Hermite 
operator B* in a proper solenoidal functional space. Of course, this belongs to classic 
theory of self-adjoint operators; see Birman-Solomjak [H p. 48]. Moreover, in a full 
capacity and specially for the NSEs in M 3 and M 2 , this theory was already developed 
(though has been used for large-time behaviour of solutions in the opposite limit t — >■ +00) 
with eigenf unctions of the adjoint (in L 2 -metric) operator B. However, since this case is 
self-adjoint, we can directly use such a theory in the complementary blow-up limit t — > 0~. 

2. Hermitian spectral theory of the linear rescaled operator B*: point 

SPECTRUM AND SOLENOIDAL HERMITE POLYNOMIALS 

Thus, approaching the point (0, 0) in the blow-up manner (11.41) . one observes Hermite's 
operator B* as the principal linear part of the rescaled equation (11.51) . Writing it in the 
corresponding divergent form, 

(2.1) B*v^V-(p*Vv), 

|y|2 

where the weight is p*(y) = e 4 > 0, we observe that the actual rescaled evolution is 
now restricted to the weighted L 2 -space L 2 , (M 3 ), with the exponentially decaying weight 
p*(y). Here, B* is the ("adjoint") Hermite operator with the point spectrum jU p. 48] 

(2.2) cr(B*) = {Afe = -f, k= |/3| = 0,1,2,.„} (/3 is a multiindex), 

where each A& has the multiplicity 2 ^ OT ^ = 3). The corresponding complete 
and closed set of eigenf unctions $* = {i^%(y)} is composed from separable Hermite poly- 
nomials. Note another important property of Hermite polynomials: 

(2.3) V^a, any derivative D 7 ip% is also an eigenfunction with k — \/3\ — I7I > 0. 



Of course, using P in (11.61) emphasizes an unpleasant fact that the NSEs are a nonlocal parabolic 
problem, so that a somehow full use of order-preserving properties of the semigroup is illusive; though 
some "remnants" of the Maximum Principle for such second-order flows may remain and actually appear 
from time to time in some results (but these are completely illusive for more difficult fluid models (|1.3I) V 
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Recall that [I] 

(2.4) polynomial set $* is complete and closed in L 2 , (IR 3 ). 

Further spectral properties are convenient to demonstrate using the linear operator B, 

(2.5) B = A + |y-V + |J in L 2 p (R 3 ), where p=j>, 

which is adjoint to B* in the dual L 2 -metric. It has the same point spectrum and the 
corresponding eigenf unctions are multiple of the same Hermite polynomials according to 
the well-known generating formula: 

(2.6) My) = t $r Dfi Hv) = rfi{y)F{y), where F{y) = J ^e-\yf 

is the rescaled kernel of the fundamental solutions of D t — A in M 3 x R + . Then, the 
bi-orthonormality holds: 

(2.7) (^,V 7 ) = Vy for any f3, 7, 

where (-, ) is the scalar product in L 2 (M 3 ). As is well known, this dual L 2 -metric can be 
also treated as a weighted one in L 2 ,(M. 3 ) (where B* becomes symmetric): 

(ty, ^ 7 ) = / R3 F( y ) ^* dy ~ (v£, 4>;) P *, 

since F(y) ~ P*{y), up to a constant multiplier. However, we prefer to keep the bi- 
orthonormality in the non-symmetric form (12. 7p . since a similar condition occurs in the 
principally non-symmetric Burnett cases; see (1B.20J) . 

Obviously, one needs to consider eigenfunction expansions in the solenoidal restriction 

(2.8) L 2 »(M 3 ) = L^(M 3 ) 3 n{divv = 0}. 

Indeed, among the polynomials $* = {ifi%}, there are many that well-suit the solenoidal 
fields. Namely, introducing the eigenspaces 

$*=Span{^, \P\ = k], k>l, 

in view of fl2.3|) . div plays a role of a "shift operator" in the sense that 



(2.9) div : $r -> n 



fe-1- 



We next define the corresponding solenoidal eigenspaces as follows: 

(2.10) SI = {v* = [v* 1 ,v;,v*] T : divv* = 0, v* e $fc}, where dim S* k = k(k + 2); 

see |221 [231 124] and further references therein. 

The study [22] deals with global asymptotics as t — > +00, where the adjoint operator 
B given in ( 12. 5 p occurs. Since B is self-adjoint in L 2 (M. 3 ), almost all the results from 
[23j Append. A] are applied to B*. For a full collection, see [51 [6] for further large-time 
asymptotic expansions and self-similar solutions. In particular, this made it possible to 
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construct therein fast decaying solutions of the NSEs on each ID stable manifolds with 
the asymptotic behaviouijj 



\ip(x,t) ~ t Xk 2 vg(-7=) + ... as t — > oo, where v^ 



v^F e S k 



^D?F(y), 



h), 



(2.11) 

are solenoidal eigenfunctions of B. Namely, taking 

(2.12) v = [v 1 ,v 2 ,v 3 ] T E S k , Vi e® k = Span {ipp - — == T - 
where F stands for the rescaled Gaussian in (12. 6p . we have that 

(2.13) div v = (vi) w + (v 2 ) y2 + {v 3 ) y3 = div (v*F) = (div v*)F - § y ■ v* F. 

This establishes a one-to-one correspondence between solenoidal eigenfunction classes SI 
in ( 12TT01 for B* and S k in (ETTIjI for B; see fl2Tnjl - (|2TTH|l below for the first eigenfunctions 
V/3 = v |(-^- Therefore, dim S k = k(k + 2), etc.; see details and rather involved proofs of 
the asymptotics ( 12. lip for fc = 1 and 2 in [22J. 

In particular, those solenoidal Hermite polynomial eigenfunctions of B* can be chosen 
as follows [23j p. 2166-69] (the choice is obviously not unique, normalization constants 
are omitted): 

Ao = 0: v* = [l,l,l] T 



(2.14) 
(2.15) 



Ai — — ^ 



'n 



e (the first solenoidal Hermite polynomial), 

(dimS* = 3); 








2/3 




-2/2 


-2/3 


, v^ 2 = 





, vj 3 = 


2/i 


Vi 




-2/i 








(2.16) 



'26 



-1 


: v 2i = 


"4 


- 2/ 2 2 - 
2/12/2 


2/1" 


, V 22 = ; 


2/12/2 
4-2/?- 


-2/1 




-2/12/3 






-2/22/3 


V24 = - 




-2/12/3 

2/12/2 


. V 25 


= — 


2/22/3 



-2/22/1 
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-2/22/3 




2/12/2 




2/ 2 2 - 2/1 




2/22/3 


> V 27 = 


2/1 - 2/? 


) V 28 = 


-2/12/2 




y\-y\ 






.- 


2/22/3 _ 






J 


'life 



'23 



2/12/3 

-2/22/3 

-2/i 2 



2/ 2 2 



(dim^ = 8), etc. 



We need the following final conclusion. By (12.41) . the set of vectors $* 3 is complete and 
closedj in L 2 , ( 



(2.17) 



VveL}. 



3 ) 3 , so that 

3\3 



v = E W c^, v;G$f, A;=|/3|>0, 



We present here only the first term of expansion; as usual in dynamical system theory, other terms in 
the case of "resonance" can contain In t- factors (q.v. pQ for a typical PDE application); this phenomenon 
was shown to exist for the NSEs in M 2 [23, p. 236]. 

Note a standard result of functional analysis: all reasonable polynomials are complete in any weighted 
L p -space with an exponentially decaying weight; see the analyticity argument in Kolmogorov-Fomin |25[ 
p. 431]. 
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where we use the vector notation 

(2.18) C/J (r) = [4(r), c%t), c%r)} T =» erf = [c%r)v^ c%r)v^ 4(r)v} 3 ] T . 
It then follows from ([2T?j) - (l2~l?]) that 

(2.19) polynomial set $* = $* 3 H {divv = 0} is complete and closed in L 2 pt (R 3 ). 

In what follows, we always assume that we deal with "solenoidal" asymptotics involving 
eigenfunctions as in ( I2.10p . 

As was promised to go in parallel, for Burnett equations in f |1.3j) . the blow-up rescaling 
and elements of linear solenoidal spectral theory are given in Appendices A and B. 

3. First application of Hermitian spectral theory: Sturmian local 
structure of zero sets of bounded solutions and unique continuation 

3.1. A dynamical system for Fourier coefficients. Consider the NSEs (11.21) . We 
assume that, in a neighbourhood of the point (x,t) = (0,0~), the solution u(x,t) is 
uniformly bounded and is such that the eigenfunction expansion of the corresponding 
rescaled function 

(3.1) u(y,r) = £g8)C/9(t)vj§(j/), where C/3 v£ = {c 1 v^ l ,c 2 v; 2 ,c 3 v} 3 ) T G L 2 p *{R 3 ), 

converges in L 2 ,(M. 3 ), and moreover, uniformly on compact subsets. These convergence 
questions of polynomial series are standard [I]; see also [HI US], where generalized Hermite 
polynomials occur and further references and details are given. In particular, if u(x,t) 
remains bounded for all t G [—1,0], then, obviously, for such bounded data uo, the 
convergence in (13.11) always takes place. 

Then, the expansion coefficients satisfy the following dynamical system (DS): 

Jc/3 = (A/?- i)c / 3 + £ (Q7) d Q7/3 c a c 7 for any \/3\ > 0, 

y where d ai p = — (P (v* • V)v*, v^) for all a, 7. 

Since (13. ip is a standard eigenfunction expansion via Hermite polynomials of a given 
bounded smooth rescaled solution u(y, r) G H 2 *(M?), the quadratic sum on the right- 
hand side converges. Recall that, moreover, according to the blow-up scaling (jl.4p . we 
actually deal with bounded and uniformly exponentially small rescaled solutions satisfying 

(3.3) |u(y,r)| < C7e"5 i n R 3 x R + . 

The DS (13. 2p is difficult for a general study. For instance, it is supposed to contain the 
answer to the existence/nonexistence of the L°°-blow-up question (The Millennium Prize 
Problem, p3])> i.e., whether there exists a Type-II blow-up at the internal point (0,0 _ ) 
(see a discussion in [201 § 5]). 

For regular points, the DS (13. 2p can provide us with a typical classification of multiple 
zeros and nodal sets of solutions. Note again that this kind of study was first performed 
by Sturm in 1836 for linear ID parabolic equations [30]; see historical and other details 
in pH Ch. 1]. 
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Thus, following these lines, we clarify local zero sets of solutions of the NSEs at regular 
points. Assume that 

(3.4) u(0,0) = 0. 

In this connection, recall that the first eigenfunction of B* with A^ = 

(3-5) v*(y) = (l,l,l) T , 

is the only ones that have an empty nodal set. Then, bearing in mind the blow-up scaling 

1 T 

term (1 — £)" = e^ in (J1.4J1 . we have to assume that 

(3.6) c o(t) = or c (t) — > as r — > +oo exponentially faster than e~5. 

3.2. Polynomial nodal sets for the Stokes equations. A first clue to a correct under- 
standing of the DS (13. 2p is given by the Stokes equations (II. ip . i.e., without the quadratic 
convection term. Then (13. 2 p becomes a linear diagonal system and is easily solved: 

(3.7) c^ = (Xp-Qcp ==► cp(r) = cp(p)e ~ for any \/3\ > 0. 

Therefore, according to (13. ip (and bearing in mind the completeness-closure of the Hermite 
polynomials), all possible multiple zero asymptotics for the Stokes problem (its local 
"micro-scale turbulence" ) are described by finite solenoidal Hermite polynomials, and the 
zero sets of rescaled velocity components also asymptotically, as r — > +oo (i.e., t — > 0~) 
obey the nodal Hermite structures. 

3.3. Nodal sets for the Navier— Stokes equations. Consider the full nonlinear dy- 
namical system ( 13. 2p . which on integration is 

(3.8) c/j(r) = ^(Oje 2 -e 2 J ^ (a7) d a7/3 (c Q c 7 )(s)e 2 ds. 



It follows that the nonlinear quadratic terms in ( I3.8p . under certain assumptions, can 
affect the rate of decay of solutions near the multiple zero. As usual in calculus, the 
indeterminacies in this integral quadratic term can be tackled by L 'Hospital rule, but 
technically this is very difficult. 

Since we are mainly interested in the study of nodal structures of solutions by using the 
eigenfunction expansion (13.11) . we naturally need to assume that it is possible to choose the 
leading decaying term (or a linear combination of terms) in this sum as r — > +00. Then 
obviously these leading terms will asymptotically describe "polynomial-like" structures of 
nodal sets as t — > 0~. For PDEs with local nonlinearities, this is done in a standard manner 
as in |16[ § 4]; in the nonlocal case, this can cause technical difficulties. However, the DS 
(13. 2p looks (but illusionary) as being obtained from a problem with local nonlinearities. 
In other words, the nonlocal nature of the NSEs is hidden in (13 ,2p in the structure of the 
quadratic sum coefficients {da^p}, and this does not affect the nodal set behaviour for 
some classes of multiple zeros. We will check this as follows: 
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Resonance zeros . We consider a "resonance class" of multiple zeros. Namely, let us 
assume there exist a multiindex subset B and a function h(r) — > such that 

c f3( r ) ~ h(r) as r — > +00 for any f3 E B, 

\cfi(r)\ <C |h(r)| as r — ► +00 for any 13 ^ B. 

In other words, only the coefficients {(^(r), /3 G B} are assumed to define the nodal set 
via (13. ip . and other terms are negligible asr-} +00. Under the natural assumption of a 
strong enough convergence of the quadratic sums in ( 13. 21) (this is expected not to be valid 
in the case of singular blow-up points only), taking the ODEs from (13. 2p for each /3 G B 
yields, for r 3> 1, 

(3.10) tfi = (A/? - |)ca + o(c^), where Cjg(r) ~ h(r). 

Hence, the asymptotic balancing of these equations must assume that, asr-y +00, 

(3.11) h~(A/j-|)h =► C/3 (r) ~ h(r) ~ eH"!)^ and |/3| = Jfc, 

where we may omit lower-order multipliers. Thus, there exists a k > 1 such that |/3| = A; 
for any /3 G £>. One can see that, for such "resonance" multiple zeros, the nonlocal 
quadratic term in (13. 2p is not important. Thus, in the resonance zero class prescribed by 
(13. 9p . as r — > +00, on compact subsets in y, similar to Stokes' problem, 



(3.12) the nodal set of u(y, r) is governed by some solenoidal Hermite polynomials. 



Polynomial structure of multiple zeros is universal . Note first that the con- 
clusion that, locally, for any zero of finite order at (0, 0), 

(3.13) nodal sets of u(x, t) are governed by finite-degree polynomials 

is trivially true for any sufficiently smooth solution. Indeed, this follows from the Taylor 
expansion of such solutions 

(3.14) u(M)= E C^x»(-ty + R K (x,t), where C^ = ^ (D^u) (0, 0) 

(M.M<*) 

and Rk = o(|x| A '(— £) ) is a higher-order remainder. Translating (I3.14p via (II. 4p into the 
expansion for u(y, r) yields some polynomial structure, so (I3.13P is obviously true. Thus, 
the principal feature of (13. 12[) is that the Hermite polynomials count only therein. 

Non-resonance zeros: a general classification . For the nonlocal problem (II .6p . 
there exist other non-resonance zeros. Indeed, let (0, 0) be a zero of u(x, t) of a finite order 
M> 1, i.e., 

(3.15) u(x,0) = E(H=M) a ^ <T ( 1 + °( 1 )) ~ x ° as x ^° (E(h=a/) \ a °\ ^°)- 
We now use the following expansion: 

(3.16) u(x, t) = u(x, 0) - u t (x, 0)(-t) + |f Utt(x, 0)(-t) 2 + ... , 
where, by (II. 6p . all the time-derivatives D^u(x,0) can be calculated: 

(3.17) u t (x, 0) = Au(x, 0) + (P(u • V)u)(x, 0) ~ x CT " 2 + (P(u • V)u)(x, 0), 
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with a natural meaning of Ax u ~ x a 2 . If the nonlocal term is negligible here and for 
other time-derivatives, i.e., 

Ut(x,0) ~x ff " 2 , u«(x,0) ~x ff " 4 ,..., 

then according to (13. 16[) this leads to a Hermitian structure of nodal sets. In fact, this 
repeats the pioneering zero-set calculations performed by Sturm (1836); see his original 
computations in [151 P- 3]. 

In general, the nonlocal term in (13.17J) is not specified by a local structure of the zero 
under consideration, so, obviously, it can essentially affect the zero evolution. For instance, 
as a hint, we can have the following zero: 

(3.18) Ut (0,0) = C^0 =» u(x ; t)~E(H=M)«^-H) = e_T (E(H=M)«^-l) ) 
where z = /_ t fi/ m ■ Hence, this nodal set is governed by the rescaled variable z, which is 
different from the standard similarity one y in (jl.4p . Of course, due to the nonlocality of 
the equation, many other types of zeros can be described. Actually, such non-resonance 
zeros can be governed by sufficiently arbitrary finite polynomials as the general expansion 
(13.141) suggests. However, there exists a countable family of "admissible" rescaled vari- 
ables. Recalling that bounded smooth solutions of the NSEs are analytic in both x and t 
(see references below), we have to have that there exists a finite K > 1 such that 

(3.19) £>fu(0,0)/0 and D s t u(0,0) = for s = 1,2, ..., K - 1. 
Therefore, close to (0, 0), the structure of such a multiple zero is given by 

(3.20) u(x,t)~E(H=M)«^ <r -H) i " = e-^(E( H =M)«^ <T -l), z = j^, 7 = § • 

Since K and M = \a\ are arbitrary positive integers, the exponent 7 in the expansion 
(13.201) can be an arbitrary positive rational number. Thus, the rescaled functions and 
variables in (I3.20J) exhaust all types of zero surfaces (points) focusing as x, t — ¥ for the 
NSEs in M?. 

Finally, the proof that zeros of infinite order are not possible for smooth non-analytic 
PDEs (and, as usual in such Carleman and Agmon-type uniqueness results, this occurs for 
u = only) is a difficult technical problem; see an example in [161 § 6.2]. For analytic in y 
solutions of the NSEs (see references and results in [81 [9j [32]) , this problem is nonexistent, 
and then in (I3.12p the degree of the solenoidal vector Hermite polynomials is always finite, 
though can be arbitrarily large. 

3.4. An application: a unique continuation theorem. Note another straightforward 
consequence of this analysis that gives the following conventional unique continuations 
result: let (13. 4p hold, (0, 0) be a resonance zeroj, and at least one component of the nodal 
set ofu(y,r) does not obey (I3.12p . Then 

(3.21) u = everywhere. 



Indeed, this is hard to check. However, for the Stokes equations (|1.1|) . as well as for any sufficiently 
smooth PDEs with local nonlinearities (see [16]), this assumption is not needed, so that such a unique 
continuation theorem makes a full sense. 
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Of course, this is not that surprising, since the result is just included in the existing and 
properly converging eigenfunction expansion (13. ip under the assumption (13. 9p . According 
to ( 13.14(1 . there exists another "funny version" of the unique continuation result: ( I3.21J) 
holds if a multiple zero is formed in a non-polynomial self-focusing of zero surfaces, or 
via a rescaled variable not available in f 1 3 . 2 (1 . but this is indeed trivial. 

For elliptic equations P(x, D)u = 0, this has the natural counterpart on strong unique 
continuation property saying that nontrivial solutions cannot have zeros of infinite order; 
a result first proved by Carleman in 1939 for P = -A + V, V e L^ c , in M 2 [7]; see [10(151] 
for further references and modern extensions. 

Thus, we are finishing a discussion of a first application of solenoidal Hermitian polyno- 
mial vector fields for regular solutions of the NSEs. We expect that, due to the DS (13. 2p . 
some "traces" of such an analysis and Hermite polynomials should be seen in the fully 
nonlinear study of u(y, r) at the singular blow-up point (0, 0), where, instead of (13. 4p . we 
have to assume that, in the sense of limsup^, 

(3.22) |u(0,0)| = +oo. 

Acknowledgements. The authors would like to thank I.V. Kamotski for his advice 
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Appendix A: Multiple zeros for Burnett equations 
A.l. Burnett equations. For both the systems (jl.3p . the blow-up scaling (jl.4p is replaced by 

(A.l) u(i,t) = (-t)-3%,r), y = j-^ TTIl 

so that the rescaled system (|1.6p takes a similar form 

(A.2) u T = H(u) = (B* - §l)u-P(u-V)u in R 3 xl + . 

The spectral theory of the given here "adjoint" operator 

(A.3) B * = -A 2 -iyV, where a{W) = {\ p = -f , |/3| = 0, 1, 2, ...} 

with eigenfunctions being generalized Hermite polynomials is available in [11] ; a solenoidal 
extension in the same lines is needed. Necessary spectral theory of the operator pair {B, B*} is 
developed below, in Appendix B. 

Therefore, under the same assumptions, the polynomial structure of nodal sets is guaran- 
teed for the corresponding Stokes-like and Burnett equations (resonance zeros) in fjl .3j) : and, 
moreover, for an arbitrary 2mth-order viscosity operator — (— A) m u therein. 

Remark: Burnett equations in a hierarchy of hydrodynamic models. The Burnett 
equations in (jl.3p appear as the second approximation (the NSEs (jl.2p being the first one) of the 
corresponding kinetic equations on the basis of Grad's method in Chapman-Enskog expansions 
for hydrodynamics. Namely, Grad's method applied to kinetic equations, by expanding the 
kernel of the integral operators involved via those with pointwise supports, yields, in addition 
to the classic operators of the Euler equations, other viscosity parts as follows: 



OG 



D t u = u t + (u • V)u = Z e 2n+1 A n {fi n Au) + ... = e^Au + e 2 /nA 2 u + ...) + ... ; 

n=0 

where e > is the Knudsen number Kn; see details in Rosenau's regularization approach, |29j . 
In a full model, truncating such series at n = leads to the Navier-Stokes equations (jl.2p (with 
Ho > 0), while n = 1 is associated with the Burnett equations in (jl.3p . 

Recall also that Burnett-type equations, with a small parameter, appeared as higher-order 
viscosity approximations of the Navier-Stokes equations, represent an effective tool for proving 
existence of their weak ("turbulent" in Leray's sense) solutions; see Lions' monograph [27, § 6, 
Ch. 1]. Note that the "Problem on blow-up/non-blow-up for Burnett equations in (|1.3p at an 
inner point" starts from dimensions N = 7; for iV < 6, there exists a unique global smooth 
L 2 -solution, |21| § 6]. It is expected that this open problem in R 7 is not easier at all than the 
classic Millennium Prize One for the NSEs in R 3 (ll.2p . In both cases, a construction (or proving 
its nonexistence) of a Type-II blow-up singularity is necessary, since, most plausibly, a Type-I 
self-similar blow-up solutions are nonexistent, |20| App. B]. 
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Appendix B: Solenoidal Hermitian spectral theory for 2?7ith-order operators 

We describe the necessary spectral properties of the linear 2mth-order differential operator 
in W N (m = 2 for the Burnett equations in fjl .3[) ) 

(B.l) B* = (-ir+ 1 A™-^yV„ 

and of its L 2 -adjoint B given by 

(B-2) B = (-ir+ 1 A- + ^yV J/ + ^/. 

As we have seen, for m = 1, (|B.1|) and (|B.2p are classic Hermite self-adjoint operators with 
completely known spectral properties, [H p. 48]. For any m > 2, both operators (jB.ip and 
(IB.2p . though looking very similar to those for m = 1, are not symmetric and do not admit a 
self-adjoint extension, so we follow [Tl] in presenting spectral theory. 

B.l. Fundamental solution, rescaled kernel, and first estimates. The fundamental so- 
lution b(x,t) of the linear poly-harmonic parabolic equation 

(B.3) u t = -{-A) m u in R N x R + 

takes the standard similarity form 

(B.4) b(x,t) = t-^F(y), y = ¥ fa. 

The rescaled kernel F is the unique radial solution of the elliptic equation with the operator 
dEl, i.e., 

(B.5) BF = -(-A) m F + 2k V • VF + it F = ° in RiV ' with I F = 1 - 

For m > 2, the rescaled kernel function -Fdyj) is oscillatory as \y\ — > 00 and satisfies |12|. 113] 

(B.6) |F(y)| < D e~ do M a in M^, where a = ^i G (1,2), 

for some positive constants D and do depending on m and N. 

B.2. Some constants. As we have seen, the rescaled kernel F(y) satisfies (|B.6|) . where do 
admits an explicit expression; see below. Such optimal exponential estimates of the funda- 
mental solutions of higher-order parabolic equations are well-known and were first obtained by 
Evgrafov-Postnikov (1970) and Tintarev (1982); see Barbatis [21 [3] for key references. 

As a crucial issue for the further boundary point regularity study, we will need a sharper, 
than given by (|B.6j) . asymptotic behaviour of the rescaled kernel F(y) as y — > +00. To get that, 
we keep four leading terms in (|B.5|) and obtain, in terms of the radial variable y 1— > \y\ > 0: 

(B.7) (-l) ra + 1 [F( 2m ) +m M f (2^i) + ,..] + -L yF > + N_ F = o fo r y > 1, 

Using standard classic WKBJ asymptotics, we substitute into (JB.7P the function 

(B.8) F{y) = y- So e aya + ... as y -»■ +00, 

exhibiting two scales. Balancing two leading terms gives the algebraic equation for a and So: 

(B.9) (_l)m (aa) 2 m -l = _L and 5q = rn{2N-l)-N > Q 
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By construction, one needs to get the root a of (JB.9P with the maximal Re a < 0. This yields 
(see e.g., 013]) 
(B.10) a= ^[- sin ( w ^) +[cos ( w JL_)]^- do + [bo (do>0 ). 

Finally, this gives the following double-scale asymptotic of the kernel: 

(B.ll) F{y) = y- So e - *^ [Ci sm(b y a ) + C 2 cos(b y a )] + ... as y = \y\ ->• +00, 

where C\^ are real constants, \C\\ + | C 2 1 ^ 0. In (jB.lip . we present the first two leading terms 
from the m-dimensional bundle of exponentially decaying asymptotics. 
In particular, for the Burnett equations in (|1.3|) in R 3 , we have 

(B.12) m = 2, N = 3: a = §, d = 3-2"T, b = 3^ ■ 2~f , and <5 = f. 

B.3. The discrete real spectrum and eigenfunctions of B. For m > 2, B is considered in 
the weighted space lA(R. N ) with the exponentially growing weight function 

(B.13) p(y) = e a ^ a > in R N , 



where a G (0, 2do) is a fixed constant. We next introduce a standard Hilbert (a weighted Sobolev) 
space of functions H 2m (R N ) with the inner product and the induced norm 

2m 2m 

(v,w) p = J p(y) Yl D^ v (y) D* w (y) dy and \\v\\ 2 p = J p{y) £ \D^v(y)\ 2 dy. 

RJV fc=0 RiV fc=0 

Then fl^QR^) C L 2 p (R N ) C L 2 ^), and B is a bounded linear operator from ^""(R^) to 
lA(M. ). Key spectral properties of the operator B are as follows [TT] : 

Lemma B.l. (i) The spectrum o/B comprises real simple eigenvalues only, 
(B.14) (7(B) = {A^ = -£, fc = |/3|=0,1,2,...}. 

(ii) T/ie eigenfunctions ipp (y) are given by 

(B.15) V/3(y) = ^^ ^F(y), /or any |/3| = fc. 

(iii) Eigenfunction subset (|B.14p is complete in L 2 (R) and m L 2 (IR). 

(iv) Tae resolvent (B — A/) -1 /or A <r(B) zs a compact integral operator in L 2 (R ). 

By Lemma IB.1| the centre and stable subspaces of B are given by 
(B.16) E c = Span{V>o = F} and E s = Span{^, |/3| > 1}. 

B.4. Polynomial eigenfunctions of the operator B*. Consider the operator (|B.1|) in the 

weighted space L 2 *(R N ), where (•, -) p * and || • \\ p * are the inner product and the norm, with the 
"adjoint" exponentially decaying weight function 

(B.17) p*( y ) = _i |J= e-«lvr>0. 

We ascribe to B* the domain iJ 2 I"(R ), which is dense in L 2 , (R ), and then 

B* : H 2 ™{R N ) ^ L 2 p ,(R N ) 
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is a bounded linear operator. B is adjoint to B* in the usual sense: denoting by (•, •) the inner 
product in the dual space L 2 (M ), we have 

(B.18) (Bv,w) = (v,B*w) for any v G H 2m (R N ) and w € Hp?(R N ). 

The eigenfunctions of B* take a particularly simple finite polynomial form and are as follows: 

Lemma B.2. (i) <r(B*) = cr(B). 

(ii) The eigenfunctions ip%(y) o/B* are generalized Hermite polynomials of degree \/3\ given by 

(B.19) ^*(y) = 1 [/ + E 0£l/2m]i ( _ A)mjy/3 l forany p _ 



VW- 



(iii) Eigenfunction subset (JB.19P is complete in L 2 , 

(iv) B* has a compact resolvent (B* — A/) -1 in L 2 *(R N ) for A <r(B* 



p 

2 
P 

(v) T/ie bi-orthonormality of the bases {i^p} and {ip*} holds in the dual L 2 -metric: 
(B.20) {ip0,^j) = <% 7 for any /3, 7. 

Remark on closure. This is an important issue for using eigenfunction expansions of solutions. 
Firstly, in the self-adjoint case m = 1, the sets of eigenfunctions are closed in the corresponding 
spaces, [1] (and we have used this in our previous NSEs study). 

Secondly, for m > 2, one needs some extra details. Namely, using (1B.20J1 . we can introduce 
the subspaces of eigenfunction expansions and begin with the operator B. We denote by L 2 the 
subspace of eigenfunction expansions v = "^cpijjp with coefficients c^ = (v,ip*) defined as the 
closure of the finite sums {Y1\b\<m cp^Pp} i n the norm of L 2 . Similarly, for the adjoint operator 
B*, we define the subspace L 2 * C L 2 *. Note that since the operators are not self-adjoint and 
the eigenfunction subsets are not orthonormal, in general, these subspaces can be different from 
L 2 and L 2 „, and particularly the equality is guaranteed in the self-adjoint case m = 1, a = j. 

Thus, for m > 2, in the above subspaces obtained via a suitable closure, we can apply standard 
eigenfunction expansion techniques as in the classic self-adjoint case m = 1. 

B.5. Solenoidal Hermite polynomials. The vector solenoidal Hermite polynomials are con- 
structed from (IB.19J) in a manner similar to that for m = 1; cf (I2.14p - (|2.16p . Namely, given a 
vector polynomial 

(B.21) v^ = [^ 1 ,^ 2 ,-,^J T ; where |ft| = |/3 2 | = ... = \fi N \ = \0\, 

it gets solenoidal provided that 

N 

(B.22) divvj=£(lfo) tt =0. 

For instance, for the Burnett case m = 2 and N = 3 (not all linearly independent ones are 
presented, normalization constants are omitted): 

(B.23) A = 0: v$ = [l,l,l] r , 

(B.24) Ai = -\ : vjj = [y 2 , -y S , 2/2] T , v* 2 = [y 3 , y 3 , - Vl } T , v* 3 = [-y 2 , yi,yi] T , 

(B.25) \ 2 = -l : v* 2l = [-y 2 l -yl,y 1 y 2 ,yiy?>] T , v* 22 = [y 1 y 2 ,-yl-yl,y2yz] T , etc. 
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(B.26) A 3 = -§ : v* 31 = [j/|, j/|,y?], v^ 2 = [yiyl,y 2 yl, -wivl + vl)i etc - 

(B.27) A 4 = -l: v: i = [^ + 4!,y 3 4 + 4!,yf + 4!] T , < 2 = [y W l y 2 yl -y 3 (yf + y|)] T , etc. 

As in the self-adjoint case m = 1, some technical efforts are necessary towards complete- 
ness/closure of generalized solenoidal Hermite polynomials in suitable spaces. We omit details. 
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